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Abstract We study interstellar dust evolution in various environments by means of chem- 
ical evolution models for galaxies of different morphological types. We start from the 
formalism developed by Dwek[15] to study dust evolution in the solar neighbourhood and 
extend it to ellipticals and dwarf irregular galaxies, showing how the evolution of the dust 
production rates and of the dust fractions depend on the galactic star formation history. 
The observed dust fractions observed in the solar neighbourhood can be reproduced by as- 
suming that dust destruction depends the condensation temperatures Tc of the elements. 
In elliptical galaxies, type la SNe are the major dust factories in the last 10 Gyr. With our 
models, we successfully reproduce the dust masses observed in local ellipticals (~ IQ^Mq) 
by means of recent FIR and SCUBA observations. We show that dust is helpful in solving 
the iron discrepancy observed in the hot gaseous halos surrounding local ellipticals. In 
dwarf irregulars, we show how a precise determination of the dust depletion pattern could 
be useful to put solid constraints on the dust condensation efficiencies. Our results will 
be helpful to study the spectral properties of dust grains in local and distant galaxies. 



1 Introduction 



The presence of interstellar dust affects strongly the spectral properies of 
galaxies, in particular in the ultra-violet (UV) band and in the far infrared 
(FIR). In the UV, the stellar light is absorbed and scattered (dust extinc- 
tion), whereas in the FIR, dust grains thermally re-emit most of the energy 
absorbed in the UV. 

However, the presence of interstellar dust can affect also the chemical com- 
position of the interstellar medium (ISM). Some chemical elements, called 
refractory (such as C, O, Fe), in the gas phase undergo dust depletion and a 
fraction of their total abundance is incorporated into dust grains. For these 
elements, the abundance measured in stellar surfaces is considerably higher 
than the gas phase abundances ([23], [45], [14]). 

In this paper, we aim at studing the effects of dust grains on the chemical 
composition of the ISM by means of chemical evolution models for galaxies 
of different morphological types. We start from the formalism developed by 
[15], and extend it to environments different than the solar neighbourhod, 
i.e. to elliptical galaxies and dwarf irregular galaxies. In particular, we focus 
on the effects that different star formation histories have on the evolution 
of the dust content of galaxies of different morphological types. 
This contribution is organized as follows. In section 2 we present the chem- 
ical evolution models. In section 3, we discuss our results for the S.N., for 
elliptical galaxies and for dwarf irregular galaxies. Finally, in Section 4 we 
present our summary. 

2 Chemical evolution of the interstellar dust 

2.1 Chemical evolution models for galaxies of different mor- 
phological types 

Chemical evolution models are used to to follow the evolution of the abun- 
dances of several chemical species and of the dust content of spirals, elliptical 
and irregular galaxies. 

In our picture, elliptical galaxies form as a result of the rapid collapse of a 
homogeneous sphere of primordial gas where intense star formation (SF) is 
taking place at the same time as the collapse proceeds. SF is assumed to 
halt as the energy of the ISM, heated by stellar winds and supernova (SN) 
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explosions, exceeds the binding energy of the gas. At this time a galactic 
wind occurs, sweeping away almost all of the residual gas. After the SF has 
stopped, the galactic wind is maintained by type la SNe, and its duration 
depends on the balance between this heating source and the gas cooling (we 
refer the reader to [33], [34]). 

In our scheme, large galaxies form the bulk of their stars and develop a galac- 
tic winds on shorter timescales than small galaxies, according to the "inverse 
wind" scenario [28]. Here, the models Lai and Hal of [34] are used for an 
elliptical galaxy of luminous mass 
respectively. 

For spiral galaxies, the adopted model is calibrated in order to reproduce 
a large set of observational constraints for the Milky Way galaxy [12]. The 
Galactic disc is approximated by several independent rings, 2 kpc wide, 
without exchange of matter between them. In our picture, spiral galaxies 
are assumed to form as a result of two main infall episodes. During the first 
episode, the halo and the thick disc are formed. During the second episode, 
a slower infall of external gas forms the thin disc with the gas accumulating 
faster in the inner than in the outer region ("inside-out" scenario, [31]. The 
process of disc formation is much longer than the halo and bulge formation, 
with time scales varying from ~ 2 Gyr in the inner disc to ~ 7 Gyr in the 
solar region and up to 20 Gyr in the outer disc (see table 1). In this paper, 
we are interested in the study of dust evolution in the solar neighbourhood 
(S.N.). For this purpose, we focus on a ring located at 8 kpc from the Galac- 
tic centre, 2 kpc wide. 

Finally, irregular galaxies are assumed to assemble from infall of protogalac- 
tic small clouds of primordial chemical composition, until masses in the 
range ~ 10^ are accumulated, and to produce stars at a lower rate than 
spirals. 

In figure 1, we show the star formation histories of the galactic models used 
in this work. 

The chemical enrichment of the ISM is due to three contributors: low and 
intermediate mass stars (i.e. with masses O.SMq <m,< 8M0); type la SNe, 
assumed to originate from exploding white dwarfs in binary systems, with 
total masses 2>Mq < minn ^ I6M0 [17], [29], and single massive stars with 
initial masses m > 8 — IOOM0 explode as core collapse SNe. 
The nucleosynthesis prescriptions arc common to all models. For massive 
stars and type la SNe, we adopt the empirical yields suggested by [16], which 
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Figure 1: Predicted SFRs as a function of time for different chemical evolu- 
tion models. Panel (a): S.N. (solid line) and a model for the outer regions of 
the Milky Way disc (dashed line); panel (b): two different elliptical galaxy 
models, i.e. model Lai (solid line) and Hal (dash-dotted line) of [34] ; panel 
(c): two different irregulars, one with continuous SF (solid line) and the 
other with 3 starbursts (dotted line). 
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are substantially based on the [44] and [22] yields, respectively, and are tuned 
to reproduce at best the abundances in the S.N. For low and intermediate 
mass stars (LIMS), we adopt the prescriptions by [41]. 
We adopt a Salpctcr [36] IMF for ellipticals and irregulars and a Scalo [37] 
IMF for the spiral model. 

It is worth to stress that the set of chemical evolution models presented in 
this paper allows us to reproduce the main observational features of local 
galaxies [8], and to account for the local metal budget [6], as well as for the 
cosmic history of metal production [9] and star formation ([4]; [7]). 

2.2 How to model dust production, dust destruction and 
dust accretion 

The chemical evolution of an element i in the dust is computed by using 
the formalism developed by [15]. Let X(iust,i{t) be the abundance by mass 
of the element i in the dust and G{t) the ISM fraction at the time t, the 
quantity Gdust,i{t) = Xdust,i ■ G{t) represents the normalized mass density 
of the clement i at the time t in the dust. The time evolution of Gfiust,i{'t) 
is calculated by means of the following equation: 
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(for a detailed explanatio of Eq. 1, see [10]. ^{t) is the star formation 
rate. The second, third, fourth, fifth and sixth terms of Eq. 1 represent the 
contributions by stellar sources of various types. These terms contain the 
quantities (^.f^, and Sj^ , which arc the condensation cfRcicncics of the 
element i in stellar winds, type la and type II SNe and regulate the stellar 
production of interstellar dust. These quantities represent the fractions of 
the element i which is condensed into dust and restored into the ISM by 
all the stellar types. For the condensation efficiencies, we use the same pre- 
scriptions as described in [15]. In general, we assume that stars can produce 
two different types of dust grains: silicate dust, composed by O, Mg, Si, S, 
Ca, Fe, and carbon dust, composed by C. The seventh and eighth terms of 
eq. 1 represent the dust destruction and accretion rates, respectively. These 
terms depend on the quantities Tdestr and Taccr, which represent the typical 
timescales for destruction and accretion, respectively. 

Dust destruction is assumed to occurr in SN shocks. Following the sugges- 
tions by [32] and [15], for a given element i the destruction timescale Tdestr 
can be expressed as: 

rdestr,i = {eMsNR)-' " ^ (2) 

RsN is the total SNe rate, including the contributions by both type la and 
type II SNe. Msnr is the mass of the interstellar gas swept up by the 
SN remnant. For this quantity, [32] suggests a typical value of Msnr ~ 
6800Mq, which is in agreement with the results from numerical studies of 
SN evolution [40]. As suggested by [32], typical values for the destruction 
efficiency e in a three-phase medium as the present-day local ISM are around 
0.2, hence we assume: 

eMsNR = 0.2 X 68OOM0 = 136OM0 . (3) 

Dust accretion is assumed to occur in dense molecular clouds, where volatile 
elements can condensate onto pre-existing grain cores, originating a volatile 
part called mantle ([15], [19]). For a given element i, the accretion timescale 
Taccr Can be expressed as: 



Taccr- = To,i/(l - /i) (4) 

where 

fi = —Q-^ (5) 
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According to eq. 4, the accretion timescale is an increasing function of the 
dust mass. For the timescale ro,i, typical values span from ^ 5 x 10^ yr, of 
the order of the lifetime of a typical molecular cloud [15], up to ~ 2 x 10^ 
yr [15]. In this paper, we assume a timescale ro.j = 5 x l(fyr. 
In elliptical galaxies, we assume that dust accretion occurs only during the 
starburst epoch, when large amounts of cold gas and molecular H are avail- 
able. 

The observed molecular H content in dwarf irregular galaxies is very small, 
with molecular-to-atomic gas fractions of ~ 10% or lower ([25], [13]). Moti- 
vated by these observational results, we assume that no accretion can occur 
in irregular galaxies. Finally, the last term of eq. 1 accounts for possible 
ejection of dust into the inter galactic medium (IGM) by means of galactic 
winds. This term is absent in the equation for the S.N. model, but is taken 
into account in the elliptical and irregular models. 



3 Results 

3.1 Dust evolution in the Solar Neighbourhood 

In Figure 2, we show the predicted evolution of the carbon (lower panel) 
and silicate (upper panel) dust production rates calculated by means of 
the chemical evolution model for the (S. N.). Concerning the C dust, its 
production is dominated by LIMS throughout most of the cosmic time. At 
the present day (Tq ~ 13 Gyr), a significant contribution is also due to type II 
SNe, with type la SNe playing a negligible role. The frequent discontinuities 
in the lines are due to the effect of the star formation threshold. 
The production of silicate dust is dominated by type II SNe throughout 
most of the time. At late times, the contributions by type la SNe and type 
II SNe are comparable. On the other hand, LIMS are negligible contributors 
to the Si dust production rate. 

In Fig. 3 we show the predicted present-day fractions in dust, i.e. the ratios 
between the amount of a given element locked into dust and its total ISM 
abundance, for the elements studied in this work, compared to the values 
observed by [23] in the Local Interstellar Cloud. It is worth to note that 
the Local Interstellar Cloud consists of a warm medium, characterized by 
a temperature of ~ 6000^, whereas our models provide the depletion of 



7 




Figure 2: Predicted dust production rates from various sources for a chem- 
ical evolution model of the S.N. In the lower (upper) panel, we show the 
results for the carbon (silicate) dust. Solid lines: contribution by low and 
intermediate mass stars (LIMS). Dotted lines: contribution by type II SNe. 
Dashed lines: contribution by type la SNe. 
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the cold gas, with temperatures lower than ~ IQOK. In this work, the 
dust fractions observed by [23] are taken as reference. Prom the logarithmic 
depletions S plotted in Figure 2 of [23], it is possible to derive the dust 
fractions of the cold gas through the formula 

/ = 1 - 10-^ (6) 

For C, O and S, the logarithmic depletions of the cold medium are very 
similar to the ones of the Local Interstellar Cloud. For Si, Mg and Fe, the 
dust fractions of the cold medium arc higher than the ones of the warm 
medium and are within the error bars plotted in Figure 3. As shown by 
[10], the dust fractions are nearly independent from the choice of the dust 
condensation efficiencies and with the assumption of a constant destruction 
efficiency e it is not possible to obtain a satisfactory fit to the observed dust 
fractions for all elements, in particular for C, O, and S. The observed dust 
fractions can be reproduced by assuming that the dust destruction efficiency 
depends on the physical properties of the chemical element, in particular on 
the condensation temperature Tg of that element. Tc express the volatility 
of the elements in dust [26]. In general, elements with higher condensation 
temperatures are more likely to aggregate into dust grains and are more 
resistant to destruction. [26] has calculated the condensation temperatures 
for various elements for a solar-system composition gas. Her results indicate 
condensation temperatures of 78 K, 182 K and 704 K for C, O and S, re- 
spectively. On the other hand, the condensation temperatures calculated for 
Fe, Si and Mg are 1357 K, 1529 and 1397, respectively, hence considerably 
higher than for C, O and S. 

As shown by Fig. 3, by assuming for C, O and S a dust destruction effi- 
ciency ec,o,S = 0.8 and for Fe, Si, and Mg eFe,Mg,Si = 0.2, it is possible to 
reproduce the observed dust fractions with good accuracy. 

In Figure 4, we show the time and metallicity evolution of the dust-to-metal 

ratio Dz (defined as the ratio between the dust mass and the metal mass 
surface density in the ISM) for the S.N. model and for for the outermost 
regions of the spiral disc. We compare our predictions to the values observed 
in the S.N. and in Damped Lyman Alpha (DLA) systems, which are believed 
to represent the high-redshift counterparts of local gas-rich galaxies, such 
as spirals and dwarf irregulars [5] . This figure shows that the model for the 
outskirts of a spiral disc presents dust-to-metal values and metallicities fully 
consistent with the values obsereved in DLAs. 



9 



Production + Destruction + Accretion 



3 

■a 



a 
o 



1 - 



o 

£ 0.5 



M, 



0.2 



Si Mg Fe 



Figure 3: Fractions in dust for various elements. The solid pentagons are the 
predicted present-day fractions calculated by assuming a dust destruction 
efficiency of e = 0.2 for all elements. The open diamonds are the present-day 
fractions calculated assuming a dust destruction efficiency of ec,o,s = 0.8 
for C, O and S and whereas for eFe,Mg,Si = 0.2 for Fe, Si and Mg. The solid 
and open squares are the fractions observed by [23] in the Local Interstellar 
Cloud using the set of cosmic abundances specified in their Tables 2 and 3, 
respectively. 
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Figure 4: Upper panel: predicted evolution of the dust-to-metal ratio in the 
S.N. (solid line) and at 16 Kpc from the Galactic centre (dashed line). The 
open square with error bar is the observational value for the Galaxy assuming 
the dust-to-gas ratio estimated by [20] and a solar metallicity of Zq = 
0.0133 [26]. Lower panel: predicted dust-to-metal ratio as a function of the 
metallicity Z in the S.N. (solid line) and at 16 kpc from the Galactic centre 
(dashed line). The dust-to-metal ratio and the metallicity are normalized 
to the Galactic value and to the solar value, respectively. The dotted region 
is the range observed in DLA systems [43]. The open square is the value 
observed in our Galaxy. 
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Figure 5: Production rates for several dust sources for C and silicates pre- 
dicted by model Hal+Itoh. 

3.2 Dust evolution in elliptical galaxies 

In this Section, we use the best values for the parameters derived above by 
studying dust in the S.N., in order to extend our formalism to the class of 
elliptical galaxies. By looking at Fig. 5, showing the predicted evolution of 
the dust production rates in an elliptical galaxy, we see that in ellipticals, 
the SNe la are the major dust factories in the last 10 Gyr. 

After the onset of the galactic wind, ellipticals are basically devoid of cold 
gas. Therefore, we assume that, during this phase, no dust accretion can 
occur in elliptical galaxies. We note that [32] estimates eMsnr ~ TOM© 
for a hot and rarefied medium as the gaseous halos surrounding ellipticals, 
therefore we run two models in which eMsNR = ISGOMq until the galactic 
wind and then we have an instantaneous transition to eMgNR ~ 70 Mq. 
The pre-wind dust evolution is obviously unaffected, whereas it leads to 
substantial changes in the late stage of the galactic evolution. According to 
the galaxy mass, we call them Lal+MK and Hal+MK, respectively. We 
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further modified the destruction treatment by implementing a thermal sput- 
tering term, which is thought to be the dominant source of dust destruction 
in hot plasmas. In particular, following [21] and [1], we assume that, in a 
~1 kcV plasma, nearly 90% of the dust grains will evaporate by thermal 
sputtering in Tciestrap,i ~ 10^/ne(yr cm"-^), where the electron density has 
been self-consistently evaluated at each timestep. This translates into a new 
destruction term, namely: 

Gdust,i ^ ^^m\/r \'^SN , G(lust,i /-,\ 
= Lrdust.ii'UiWQj 1 {1} 

'^destr (^gas 'TdestrsP,i 

According to the galaxy mass, we call the models featuring this particular 
term as Lal+Itoh and Hal+Itoh, respectively. When studying chemical evo- 
lution of the dust in ellipticals, it may be important to investigate whether 
the inclusion of the dust may help in solving the so-called iron discrepancy 
[1]. The ASCA satellite provided the first reliable measures of the iron abun- 
dance in the hot ISM of ellipticals (e.g. [2]; [27]), finding values much lower 
than solar, at odds not only with theoretical models for elliptical galaxies 
available at that time ([1]; [30]), which predicted that their ISM should ex- 
hibit [Fe/H]> 0, but also with the mean metallicity of the stellar component 
inferred from optical spectra. By considering the most recent data, the Fc 
discrepancy is partially alleviated, but clearly still present when comparing 
the chemical evolution results by [34], predicting [Fe/H]> 0.85 in the ISM 
of elliptical galaxies, to the abundance measurements in the X-ray spectra, 
indicating emission-weighted Fe abundances up to ~ 2 — 3 solar in a sample 
of 28 early-type galaxies, corresponding to [Fc/H]=0.3 - 0.48 [18]. 
Fe condensation in dust is among the possible physical mechanisms often 
invoked to solve this issue (see [1] for a comprehensive analysis). Recent far- 
infrared observations, in fact, claim that the dust mass in ellipticals could 
be ~ 10^~ ''Mq, a factor of ten higher than previous estimates [39], with 
~ 2 • 1O^M0 of dust residing in the very central galactic regions [24]. In- 
dependent recent observations carried out by means of the SCUBA camera 
on local ellipticals [42] seem to confirm these values, although the sample 
should be enlarged and the effects of synchrotron radiation clarified. These 
new observational values are well reproduced by our models, which predict 
a dust mass of ~ 1.1 • 10^ for more than 10 Gyr old giant ellipticals 
{rsj 0.8 • 10^ Mq of dust for a more typical galaxy like our model Lai). 
In Figure 6, we show the predicted evolution of the Fe abundance in the 
hot phase of a giant elliptical galaxy by adopting the [21] (solid line) and 
the [32] (dotted line) prescriptions. The best results have been obtained by 
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Figure 6: Evolution of the Fe abundance in the hot phase of an ellipti- 
cal galaxy for the Lai model. Solid lines: Lal+Itoh model; Dotted Line: 
Lal+MK model. 

adoptiong the Mckee presciptions for dust destruction. With these prescrip- 
tions, the predicted prcscnt-day iron abundance in the hot gas in the case 
of a and a giant elliptical is [Fe/H]=+0.3 and +0.15 dex, respectively. 
These values are in excellent agreement with the available observations. 

Moreover, since the Fe emission line dominates the X-ray spectrum of the 
hot ISM, a reduced Fe abundance in the gas leads to a better agreement be- 
tween the predicted X-ray luminosity (2.1 • lO^^crg scc^^) and temperature 
{r^ 0.7 kev) for the hot halo of giant ellipticals with the observed values, 
than in [34] case. 

To conclude, we suggest that the iron discrepancy might be solved by tak- 
ing into account the dust condensation of Fc. New and elliptical-dedicated 
calculations of the sputtering by hot (i.e. T > lO^K) plasmas are needed to 
eventually asses this issue. 
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3.3 Dust evolution in Irregular galcLxies 



In Figure 7, we show the predicted dust production rates as functions of 
time for the models of an irregular galaxy with continuous SF (IC) and a 
starburst irregular. 

All the dust production rates calculated for the IC model have a smooth be- 
haviour and in general present a continuous increase, lasting several Gyrs, 
and then reach a plateau. On the other hand, the production rates calcu- 
lated for the starburst irregular reflect the shape of the SF and present a 
gasping or intermittent behaviour. An interesting difference between the 
dust production rates predicted for irregular galaxies and the one predicted 
for the S.N. model concerns the production of silicate dust by low and inter- 
mediate mass stars. In the case of the S.N. model, the Si dust production 
by LIMS starts immediately after the first SF episode. On the other hand, 
for the IC model Si dust production by LIMS starts at ~ 5 Gyr and for 
the starburst model it does not start at all. This depends on the adopted 
prescriptions for the dust condensation efficiencies and on the C/0 ratio in 
dwarf irregulars, as can be seen in [10]. 

When studying dust evolution in irregular galaxies, it may be interesting 
to study the dust-to-gas and dust-to-metal ratios in order to compare the 
predictions with some directly-measurable quantities. The dust-to-gas ratio 
is defined as: 

where M^ust and Mhj are the dust and HI mass determined from the far 
infrared and from the radio emission, respectively [25]. In the models, the 
neutral H mass is Mhi = Mgas • Xh, where Mgas is the total gas mass and 
Xh is the H mass fraction. In Figure 8, we show the D' calculated for the IC 
model and for the starburst model as a function of time (upper left panel) 
and as a function of the metallicity (lower left panel), and compared to a 
set of observations in local dwarf galaxies. At any time, the IC model is 
characterised by higher D' values than the starburst model. 
In the lower panels of Figure 8, as a proxy of the metallicity we use the 
quantity 12-|-log(0/H). The observational values used here are taken from a 
compilation by [25] (sec caption of Fig. 8), and refer to BCD (solid circles) 
and dirr (open squares) galaxies. We note that the predicted present-day 
values are in general higher than the majority of the observed values. The 
compilation of observations presented by [25] is based on IRAS data. The 
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Figure 7: Predicted dust production rates from various sources for a chem- 
ical evolution model of an irregular galaxy with continuous SF (left) and a 
starburst irregular galaxy (right). In the lower (upper) panels, we show the 
results for the carbon (silicate) dust. Solid lines: contribution by low and 
intermediate mass stars (LIMS). Dotted lines: contribution by type II SNe. 
Dashed lines: contribution by type la SNe. 
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observed determinations of the dust mass are likely to represent lower limits 
to the actual values, owing to the undetected cold dust component ([15], 
[35]). [25] estimate a total error in the dust-to-gas ratio to be of a factor of 
4, taking into account all the possible factors of uncertainty (contribution 
from very small grains, cold dust, molecular gas and variations in the HI and 
optical diameter). By taking into account this uncertainty (solid and open 
inverted triangles in Fig. 8), the observed dust-to-gas ratios are consistent 
with our predictions for both irregular types. 

The dust-to-metals ratio D'^ observed in dwarf irregulars can be defined as: 



^ Mhi -Z Z ^ ' 

In dwarf irregular galaxies, dust accretion is assumed to be negligible, hence 
the dust fractions are determined from the balance between the dust con- 
densation efficiencies and the dust destruction process. In the right upper 
and lower panels of Figure 8, we show the predicted evolution of I?^ in dwarf 
irregular and starburst galaxies as a function of time and mctallicity, respec- 
tively, calculated assuming the set of condensation efficiencies as suggested 
by [15] (thick lines) and a constant value of 5f^ = Sj"' = = 0.1 (thin 
lines). The adoption of different condensation efficiencies has noticeable ef- 
fects on the dust-to-metal ratio. In particular, the D'^ values calculated 
by assuming the condensation efficiencies by [15] are larger than the ones 
predicted by assuming a constant value of 0.1. 



In the lower right panel of Fig. 8, we show the predicted evolution of D'^ as 
a function of the mctallicity for the IC model (solid lines) and for the star- 
burst model (dashed lines). The big open circles and the big open squares 
are the predicted present-day values for D'^ for the starburst and IC model, 
respectively. In the S.N., the dust fractions are determined by the balance 
between dust destruction and dust accretion, hence they are independent 
on the dust condensation efficiencies. Fig. 8 shows that, in principle, the 
measure of the dust-to-metal ratio or of the dust depletion pattern in dwarf 
irregulars could be helpful to put solid constraints on the dust condensation 
efficiencies. 
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Figure 8: Left: Evolution of the dust-to-gas ratio calculated for the IC 
model (solid lines) and for the starburst model (dashed lines) as a function 
of time (upper panel) and metallicity (lower panel). In the lower panel, the 
solid circles and open squares are D' values observed in local dirr and BCD 
galaxies, respectively [25]. The solid and empty inverted triangles represent 
upper limits to the observations for local dIrr and BCD galaxies, respectively, 
calculated by assuming an uncertainty of a factor of 4. Right: Evolution 
of the dust-to-metal ratio calculated for the IC model (solid lines) and for 
the starburst model (dashed lines) as a function of time (upper panel) and 
metallicity (lower panel). The thick and thin lines represent the predicted 
D' calculated assuming for the condensation efficiencies the values suggested 
by [15] and a constant value of 0.1, respectively. All the symbols are as in 
the left lower panel. 
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4 Summary and future perspectives 



By means of chemical evolution models for galaxies of different morpholog- 
ical types, we have performed a study of the cosmic evolution of the dust 
properties in different environments. By starting from the same formal- 
ism as developed by [15], we have carried on a deep study of the space of 
the parameters used to model dust evolution and, thanks to the uptodate 
observations available in the solar vicinity, performed a fine tuning of the 
parameters. We have extended our study to ellipticals and dwarf irregular 
galaxies for which, for the first time, dust evolution has been calculated 
by means of chemical evolution models relaxing the instantaneous recycling 
approximation, namely taking into account the stellar lifetimes. Our main 
results can be summarised as follows. 

1) The evolution of the dust production rates and of the dust fractions are 
strongly determined by the SF history. 

2) We can reproduce the dust fractions observed in the solar neighbourhood 

by assuming that the dust destruction efficiency depends the condensation 
temperatures Tc of the elements. In general, more volatile elements have 
lower Tc and higher dust destruction efficiencies. 

3) In elliptical galaxies, type la SNe are the major dust factories in the last 
10 Gyr. With our models, we successfully reproduce the dust masses ob- 
served in local ellipticals (~ 10^ M©) by means of recent FIR and SCUBA 
observations. 

4) We have shown that, by considering a reduced dust destruction in a hot 
and rarefied medium, dust is helpful in solving the iron discrepancy observed 
in the hot gaseous halos surrounding local ellipticals [1]. In this medium, 
we predict a Fe abundance of [Fe/H]=0.15-0.3, whereas without dust [34] 
found [Fe/H]> 0.85. 

5) The two models used to study dust evolution in dwarf irregular galaxies 
present very different dust production rates. For the IC model, the pro- 
duction rates have a smooth behaviour throughout the whole evolution of 
the system. For the starburst model, the evolution of the production rates 

reflects its intermittent SF history. 

6) The predicted present-day dust-to-gas ratios for irregular galaxies are in 
general higher than the majority of the observed values. However, the de- 
termination of the dust mass in these galaxies is affected by several sources 

of uncertainty. By assuming a factor of 4 uncertainty in the observed data, 
as suggested by [25] , the observed dust-to-gas ratios are consistent with our 
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predictions. 

7) In dwarf irregulars, dust accretion is likely to play a negligible role or to 
be absent, hence the variation of the condensation efficiencies has important 
effects on the dust-to-metal ratio. In principle, a precise determination of 
the dust-to-metal ratios or of the dust depletion pattern in dwarf irregular 
galaxies could be helpful to put solid constraints on the dust condensation 
efficiencies. 

In the future, our results will be very useful for the study of the spectral 
properties of dust grains in galaxies of different morphological types. By 
combining the chemical evolution results described here with a spectropho- 
tometric code that includes reprocessing by graphite and silicate dust grains 
(GRASIL, [38]), it will be possible to analyse the spectral energy distribu- 
tions of galaxies of different morphological types (Schurer et al. 2008, in 
prep.) and to considerably improve our understanding of interstellar dust. 
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